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Abstract 


Many methods have been developed to aid analysts in identifying component modes which contribute significantly 
to component responses. These modes, typically targeted for dynamic model correlation via a modal survey, are 
known as target modes. Most methods used to identify target modes are based on component global dynamic 
behavior. It is sometimes unclear if these methods identify all modes contributing to responses important to the 
analyst. These responses are usually those in areas of hardware design concerns. One method used to check the 
completeness of target mode sets and identify modes contributing significantly to important component responses 
is mode participation. With this method, the participation of component modes in dynamic responses is quantified. 
Those modes which have high participation are likely modal survey target modes. Mode participation is most 
beneficial when it is used with responses from analyses simulating actual flight events. For spacecraft, these 
responses are generated via a structui^ dynamic coupled loads analysis. Using MSC/NASTRAN DMAP, a method 
has been developed for calculating mode participation based on transient coupled loads analysis results. The 
algorithm has been implemented to be compatible with an existing coupled loads methodology and has been used 
successfully to develop a set of modal survey target modes. 
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Nomenclature 


Abbreviations 


Matrices 


Set Notation 


DOF Degree-of-freedom 
DMAP Direct Matrix Abstraction 
Program 

DRM Data Recovery Matrix 
LeRC Lewis Research Center 
MEM Modal Effective Mass 
MPF Mode Participation 
Fraction 

NASA National Aeronautics and 
Space Administration 
WP-02 Work Package Two 
WP-04 Work Package Four 


B Damping 

k,K Stifftiess 

m,M Mass 

p,P Applied loads 

R Recovered responses 

T Transformation 

u,U Physical displacements 

u,U Physical accelerations 

<t>,d> Mode shapes 

^ Modal displacements 

• 

I Modal velocities 

\ Modal accelerations 


a a-set (analysis DOF) 

A Acceleration responses 

D Displacement responses 

e Elastic modal DOF 

h h-set (system modal DOF) 

q q-set (generalized DOF) 

r Rigid-body modal DOF 

t t-set (total physical 

interface DOF) 

T Total responses 


Introduction 


Many methods have been developed to aid analysts in identifying component modes which contribute significantly 
to component responses. These modes, typically targeted for dynamic model correlation via a modal survey, are 
known as target modes. Most methods used to identify target modes are bjised on component global dynamic 
behavior. Such methods include modal effective mass [1] and effective interfa<» mass [2], Because they are based 
on general dynanuc characteristics of a component, sometimes it is not clear if these methods identify all modes 
which contribute significantly to responses important to the analyst. Many times, these responses include those which 
are in areas of hardware design concerns. One method used to check for the completeness of target mode sets and 
identify modes contributing sig nifican tly to important component responses is mode participation. With this method, 
the participation of component modes in dynamic responses is quantified. Those modes which have high 
participation contribute significantly to responses and are likely modal survey target modes. Mode participation is 
most beneficial when it is used with responses from analyses simulating actual flight events. For components such 
as Space Shuttle payloads, these responses are gener^ed via a structural dynamic coupled loads analysis. 
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Structural dynamic coupled loads analyses for Space Shuttle/payload systems have been performed at NASA LeRC 
for many years [3]. Such analyses have supported Space Station hardware design efforts including the WP-02AVP-04 
Combined Cargo Element shown in Fig. 1. To test-vaify the dynamic model representing this Space Shuttle 
payload, a modal survey was planned. As part of the modal survey pre-test analysis activities, target modes were 
selected. The initial criteria for target mode selection was based on total modal effective mass. Experience has 
shown that this method is appropriate when the goal of a modal survey/model correlation activity is to obtain an 
analytical model which accurately predicts Space Shuttle/payload interface loads during transient analyses. For the 
Combined Cargo Element, attention was also focussed on areas of the hardware design where margins-of-safety were 
low. In particular, these areas included interfaces between major hardware components. Given the set of target 
modes selected using total modal effective mass, analysts questioned whether or not all modes contributing 
significantly to important component responses had been included. Analysts wanted payload modes contributing 
significantly to responses in areas of concern to be selected as target modes to better ensure accurate predictions of 
dynamic loads. To help analysts determine component modes important to responses in areas of concern, mode 
participation analysis was identified as an appropriate tool. The original set of target modes identified using modal 
effective mass would be modified to include component modes deemed important based on mode participation 
analyses. 

To perform mode participation analyses for Space Shuttle payloads in support of modal survey activities, a 
MSC/NASTRAN DMAP sequence was developed. The DMAP is used to quantify component mode participation 
based on transient coupled loads analysis results. Developed for multi-level MSC/NASTRAN superelement payload 
models such as that representing the WP-02/WP-04 Combined Cargo Element, a major design feature of the new 
analytical tool is that it allows for relating superelement responses to the modes of an assembled payload model. 
This is accomplished through the use of component model response DRMs. The method has been implemented to 
be compatible with the NASA LeRC coupled loads methodology and has been used successfully to develop a set 
of modal survey target modes for the WP-02/WP-04 Combined Cargo Element. 

The objective of this work was to develop a methodology that would quantify the participation of component model 
normal modes in coupled loads transient responses. In the following section, the general theory for mode 
participation calculations using component response DRMs is presented. The implementation of the general theory 
w ithin a MSC/NASTRAN DMAP sequence is detailed in a subsequent section. Finally, the sqjplication of the new 
mode participation DMAP sequence in selecting a set of Space Station WP-02/WP-04 Combined Cargo Element 
modal survey target modes is described. 


Theory 

To identify component modes which contribute significantly to transient analysis component responses, a mode 
participation methodology has been developed which uses component model DRMs. Component DRMs transform 
the responses of component model interface and generalized DOF to general component responses (accelerations, 
forces, stresses, etc.). To better understand the mode participation methodology, it is worthwhile to review 
calculating component model interface and generalized DOF transient responses and performing transformations made 
possible by using DRMs. 

Let a system for which a structural dynamic coupled loads analysis is to be performed be an assembly of many 
component models. Neglecting component damping, the equations of motion for the analysis (a-set) DOF of one 
of the component models are 


[mj{uj + [kj{uj = {pj 


( 1 ) 


where [mj, [kj, and {pj are the dynamically reduced component a-set DOF mass, stiffness, and applied loads 
matrices, respectively. The component a-set DOF responses are the displacements {u,} and its derivatives. For a 
Craig-Bampton [4] reduced component model, the a-set DOF are comprised of the interface (t-set) DOF and 
generalized (q-set) DOF. The q-set DOF correspond to the component model normal modes of vibration when all 
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t-set DOF are fixed. Given the reduced a-set size matrices of all component models, the a-set size matrices of all 
models are assembled to form the coupled system a-set DOF equations of motion: 


[MJ{UJ + [KJ{UJ = {P,} 

In Eq. (2), coupled system damping has yet to be defined, and [MJ, [KJ, and {P^} are the coupled system mass, 
stiffness, and applied loads matrices, respectively. The coupled system a-set DOF responses are the displacements 
{Uj} and its derivatives. 

Typically during a coupled loads analysis, the coupled system equations of motion shown by Eq. (2) are transformed 
to modal space. The reduction uses the coupled system free-free modes [‘b J generated via a system-level 
eigenvalue analysis to transform the n^ number of a-set DOF to n^ number of modal (h-set) DOF. Since the number 
of h-set DOF is usually much less than the number of a-set DOF (n^ « n,), the reduction is made to take advantage 
of efficiencies gained when solving for the coupled system responses using modal analysis methods. For a modal 
transient analysis, the responses for the n^ number of coupled system modal DOF are found from solving 

mjil) + [KJ{^} = {P,} 

where [MJ, [BJ, [KJ, and {PJ are the coupled system modal mass, damping, stiffness, and applied loads 
mafrirpR , respectively, and and its derivatives are the coupled system modal DOF responses. Damping matrix 
[Bh,] is generated external to the coupled system reduction of Eq. (2). Typically, [Bi^,] is a diagonal matrix with 
terms proportional to the coupled system eigenvalues. Because the system r^uction is made on all component a-set 
DOF, the coupled system modal DOF responses are functions of all components which are part of the coupled 
system. 

After solving for the coupled system modal DOF responses via Eq. (3), the a-set DOF responses for any component 
which is part of the system can be solved for. Component a-set DOF accelerations and displacements can be 
recovered using two basic methods: mode displacement or mode acceleration. For the mode displacement method 
of recovery, the component a-set DOF responses are calculated as 


{uj = 

and 

{Ua) = 


(4) 


(5) 


where [<j>J is the partition of the coupled system mode shapes [OJ corresponding to the component a-set DOF. 
For the mode acceleration method, the component a-set DOF responses are 


{u»} = 


( 6 ) 


and 


{uj = + [ua'KPj 

where "r" and "e" refer to the n^ and n, number of coupled system rigid-body and elastic modal DOF, respectively. 
In Eq. (7) [3], [<|) J and [<|) J are the rigid-body and elastic partitions of [<j) J, is a diagonal matrix containing 
the inverse squares of the coupled system elastic circular natural frequencies, [u,‘] are static displacements of the 
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component a-set DOF due to unit static physical loads applied to the coupled system, and {Pa} are tlie loads applied 
to the coupled system. Whether the component a-set DOF responses are recovered using a mode displacement or 
mode acceleration method, they are functions of the coupled system modal DOF responses. 

As explained previously, the component a-set DOF can be partitioned into n, number of total physical interface (t-set) 
DOF and n, number of generalized (q-set) DOF. Partitioning the component a-set DOF accelerations and 
displacements. 




(8a,b) 


Given the partition of the component a-set DOF and considering the partition of the coupled system modal DOF, 
the mode displacement solution for the component a-set DOF accelerations and displacements shown by Eq. (4) and 
Eq. (5), respectively, can be rewritten as 


and 





[{u,} 




(9) 


( 10 ) 


Similarly, the mode acceleration solution for the component a-set DOF accelerations and displacements shown by 
Eq. (6) can be rewritten as Eq. (9), and the displacements shown by Eq. (7) can be rewritten as 



[<I>J 




[Ut‘] 

[U^ 


a,) 


( 11 ) 


The above discussion focusses on the solution of component model a-set DOF accelerations and displacements using 
either a mode displacement or mode acceleration method. Given these responses, the efficient recovery of general 
component responses using DRMs is next considered. A DRM can transform the component a-set DOF responses 
to general component responses as 


|{Ra} 


.rriJ rr^lj [{VJ 


( 12 ) 


where (Ra) and (Rd) are the n^ number of recovered accelerations and number of recovered displacement- 
dependent responses (displacements, forces, etc.), respectively. Note that Eq. (12) shows a DRM for a component 
mode acceleration transformation to generate internal responses. Matrix [TAa] is the partition of the DRM 
transforming component a-set DOF accelerations into internal accelerations. Matrices [T^J and [Tot] combine to 
transform the component a-set DOF accelerations and t-set DOF displacements into internal displacement-dependent 
responses. Partitioning the component a-set DOF terms of Eq. (12) into t-set DOF and q-set DOF terms, Eq. (12) 
can be rewritten as 
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( 13 ) 


where 


|{Ra}| 

{R1} + {Ra} 


riAti iTa^ PatI 


{R^} + {R’} 


it'd.] [Tj,] 


{ii,} 

{%)■ 


iK) 


rrj [OaJ [oj 

iK) 


[Oa.1 [Oa.1 

iK) 


it'd.] VW ITd,] 



[Od.] ITd,] [Od.] 


(uj 

- {«,} 

HI 


(14) 


The above discussion centers on calculating component a-set DOF responses and then calculating general component 
responses using DRMs. Mode participation analysis is used to identify which component modes contribute 
significantly to component responses. The component modes are those used to reduce the component model and 
correspond to the component q-set DOF. For this application, mode participation analysis identifies important modes 
by calculating mode participation fractions (MPFs). Mode participation fractions are defined as follows: Given a 
component response at an instant in time, a particular component mode MPF is defined as the response due to the 
corresponding component q-set DOF divided by the sum of all component q-set DOF responses. In this way, 
component modes which contribute most to a response are identified by q-set IX)F with high MPFs. 

To develop an equation with which MPFs can be calculated, first consider the icomponent q-set DOF partitions of 
Eq. (14): 





{Rg} 


[TdJ. 


Equation (15) can be written in a simpler form as 

{R,} = [TgCti,} 


(15) 


(16) 


where {Rj} contains the sums of all component q-set EXDF contributions to the nj number of total recovered 
responses (nf=nA+no), and [TtJ is the partition of the DRM transforming the q-set DOF accelerations to the q-set 
DOF contributions. For any jth response within {Rf}, the sum of all q-set DOF contributions to that response is 

Rj = [Tg{5,} 


where [Tj,] is the jth row of Uxq]. The multiplication shown by Eq. (17) can also be written as the summation 


r, = er] = et;k 




(18) 


The n, number of terms shown in Eq. (18) are the individual q-set DOF responses which contribute to the jth 
response. Given these individual contributions, the MPF for the ith q-set DOF is calculated as 
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(19) 


MPF 


i _ r; _ 


R. 


R. 


As can be seen from the above derivation, calculating MPFs for a particular component response is relatively straight 
forward. The MPFs are used to identify which component modes contribute most to responses of interest. It is 
important to note that identifying modal survey target modes based upon mode participation calculations requires 
engineering judgement. Judgement is required when selecting which responses to investigate and interpreting 
anal ytical results. An illustration of this is given via an example in a following section. 


Im plementation using MSC/NASTRAN DMAP 

A mode participation methodology has been implemented as a MSC/NASTRAN DMAP sequence compatible with 
the NASA LeRC coupled loads methodology. A DMAP listing is provided in the Appendix. The DMAP sequence 
is operated as a data recovery restart using modal transient analysis results. The input data to the DMAP sequence 
are illustrated in Fig. 2. System modal DOF transient responses, system mode shapes, component model data and 
transformation matrices, and component response DRMs are required. The loads applied to the coupled system are 
also needed if system modal DOF mode acceleration recovoy is desired. Although the DMAP sequence has been 
designed to be used with the NASA LeRC coupled loads methodology, the same input data could be obtained from 
standard MSC/NASTRAN solution sequences. 

The mode participation DMAP sequence was written under the assumption that the component response DRM used 
for the calculations is of a form shown in Eq. (12). The DRM is for the recovery of internal responses using 
component model mode acceleration techniques. This DRM format is one used often within the NASA LeRC 
coupled loads methodology. However, minimum modifications could be made to the mode participation DMAP to 
accommodate different DRM types. A flow-chart illustrating the basic mode participation DMAP operations is 
shown in Fig. 3. Given system modal DOF transient responses found via Eq. (3) and component model data and 
transformation matrices, the component a-set DOF accelerations and t-set DOF displacements are calculated using 
either mode displacement (Eq. (9) and Eq. (10)) or mode acceleration (Eq. (9) and Eq. (11)) data recovery. The 
method of recovery is an option to the analyst. The component response DRM is then read in INPUTT4 format 
The DRM is created using custom DMAP sequences [3]. For a multi-level MSC/NASTIL\N superelement model, 
the DRM is used to relate component responses within any superelement to the a-set DOF responses of the 
assembled model. In this way, any general component respcmse of a multi-level superelement model can be related, 
via the DRM, to the q-set DOF responses, hence the component modes, of the entire assembled model. 

To perform mode participation analysis for one or all component response(s) defined within the DRM is an option 
to the analyst. Given that mode participation calculations can be performed for any output time instant of a transient 
analysis, the time instant for mode participation analysis is also an option. The analyst can allow calculations to be 
automatically performed for the time instant of absolute maximum response, or the analyst can input the time instant 
of interest using a parameter. Once the component response(s) and analysis time instant have been defined, the 
coixesponding row(s) of the DRM and column of the component boundaty responses are extracted. Intermediate 
calculations include individual component q-set DOF response contributions and the sum of all component q-set DOF 
response contributions shown by Eq. (18). Lastly, MPFs for the particular response(s) are generated using Eq. (19). 
These data are written to the MSC/NASTRAN .F06 file for post-processing. 


Application to Modal Survey Target Mode Selection 

A mode participation analysis methodology for general component responses was developed in support of Space 
Station WP-02AVP-04 Combined Cargo Element modal survey pre-test analysis activities. The Combined Cargo 
Element was designed to be lifted into space by the Space Shuttle. A design for the payload is shown in Fig. 1. 
In the figure, the payload is appropriately aligned with the Space Shuttle X„ Y„ and Z„ axes. The payload was 
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approximately forty-five feet long (X„ direction), fifteen feet wide (Y„ direction), and fifteen feet high (Z„ direction). 

It was designed to mount within the Space Shuttle cargo bay at six attach points (eight DOF). The payload was to 
attach in two DOF in the X„ direction, two DOF in the Y^, direction, and four DOF in the Z, direction. The finite 
element model representing the payload for transient dynamic analyses is shown in Fig. 4. The model had a weight 
of 35,400 lb and consisted of approximately 79,000 MSC/NASTRAN g-set DOF. For dynamic analyses, the model 
was substructured into a multi-level superelement model. The payload model superelement tree is shown in Fig. 5. 
Fixed at the eight Space Shuttle interface DOF, the first fifty-six payload frequencies and their translational modal 
effective masses were as listed in Table 1. 

The criteria originally used to identify Combined Cargo Element modal survey target modes was based solely on 
total translational modal effective mass. The goals for total modal effective mass were equal to 90% in the Xo 
direction, 90% in the Y„ direction, and 75% in the Z<, direction. The lower value for the Z^ direction was due to the 
redundancy in Zo attach DOF to the Space Shuttle. Choosing translational modal effective mass percentages for each 
payload fixed-interface mode listed in Table 1 greater than 2.0%, the original set of modal survey target modes was 
as listed in Table 2. Twelve payload modes were chosen. Only the modal effective masses greater than 2.0% for 
each mode are shown in Table 2. 

Modal effective mass was appropriate for identifying target modes important for accurately predicting Space 
Shuttle/payload interface loads. However, analysts also wanted the target mode set to include modes that contributed 
significantly to dynamic transient responses in areas of low payload hardware design margins-of-safety. Mode 
participation analyses were performed for key payload responses to investigate whether or not additional target modes 
could be identified. A typical area of concern was as highlighted in Fig. 4. Regions such as this corresponded to 
major hardware interfaces and were represented in the loads analysis model by the interfaces between superelements 
1 and 4 and superelements 2 and 4 (refer to Fig. (5)). Since the regions of interest were very similar in terms of 
design and dynamic response, the remaining discussion focusses on the interface between superelements 1 and 4. 
For mode participation analysis, a DRM was generated for recovering the interlace loads between superelements 1 
and 4 in the area of concern. For interface load recovery, the DRM was post-multiplied by the a-set DOF responses 
of the assembled payload model. The q-set DOF respcmses corresponded to the payload fixed-interface modes, some 
of which are listed in Table 1. As an example of the mode participation ainalyses, one of the interface loads 
investigated was aligned along the Space Shuttle Y^ axis. Significant payload q-set DOF contributions to that load 
are shown in Fig. 6. Referring to Fig. 6, payload modes 28, 29, and 30 contributed significantly to the response. 
This observation was also noted for mode participation analyses made for other interface loads. Hence, these three 
payload modes were added to the set of modal survey target modes as shown in Table 2. Note that the modal 
effective masses for modes 28, 29, and 30 are not significant with respect to the 2.0% criteria. However, mode 
participation analyses showed that these modes contributed significanfly to dynamic responses in areas of concern 
for hardware design margins-of-safety. 

In most engineering applications, a significant number of component modes may be identified via mode participation 
analyses as being significant contributors to component responses. This set of modes will probably be much larger 
than a set of target modes selected via modal effective mass or other techniques. Many times, constraints placed 
on analysts and modal survey test personnel require that the total set of target modes be prioritized. Hence, it is 
important that engineering judgement play an important role in selecting additional target modes based on mode 
participation analyses. 


Summary 

A mode participation analysis methodology has been developed and implemented using MSCl/NASTRAN DMAP. 
The methodology, designed for transient dynamic coupled loads analyses, quantifies which fixed-interface component 
modes participate most in component responses. Knowing which component modes significantly participate in 
internal component dynamic responses can aide the analyst during modal survey pre-test analysis activities in 
selecting canriidate target modes. The methodology takes advantage of efficiencies gained using MSC/NASTRAN 
superelement capabilities and component data recovery matrices. The methodology is designed to be compatible with 
the NASA LeRC coupled loads analysis methodology. 
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Table 1. WP-02AVP-04 Combined Cargo Element Fixed-interface Normal Modes Summary 
Natural Frequencies and Modal Effective Mass 


Mode 

Freq. 

X„MEM 

YoMEM 

Z„MEM 

Mode 

Freq. 

X„MEM 

Y<,MEM 

Z„MEM 


(Hz) 

(%) 

(%) 

(%) 


(Hz) 

(%) 

(%) 

(%) 

1 

6.66 


87.5 

0.1 

29 

19.09 

0.0 

0.0 

0.5 

2 

8.30 

49.8 

0.3 

21.4 

30 

19.16 

0.0 

0.0 

1.4 

3 

9.40 

46.4 

0.0 

26.7 

31 

19.41 

0.0 

0.0 

0.0 

4 

10.60 

0.2 

6.0 

0.1 

32 

19.46 

0.0 

0.0 

0.1 

5 

12.66 

0.2 

1.4 

0.3 

33 

19.51 

0.0 

0.0 

0.0 

6 

13.08 

0.0 

0.0 

2.1 

34 

19.81 

0.0 

0.0 

0.0 

7 

13.46 

0.1 

0.6 

0.2 

35 

19.86 

0.0 

0.1 

0.1 

8 

13.93 

0.0 

0.0 

0.5 

36 

19.98 

0.0 

0.5 

0.3 

9 

14.08 

0.0 

0.2 

0.0 

37 

20.06 

0.0 

0.0 

0.0 

10 

14.17 

0.0 

0.0 

0.2 

38 

20.50 

0.0 

0.0 

0.0 

11 

14.18 

0.1 

0.0 

0.2 

39 

20.82 

0.1 

0.0 

0.0 

12 

14.26 

0.0 

0.0 

4.7 

40 

20.85 

0.0 

0.0 

0.0 

13 

14.33 

0.6 

0.1 

8.6 

41 

20.86 

0.1 

0.0 

0.0 

14 

14.51 

0.0 

0.0 

3.7 

42 

20.92 

0.0 

0.0 

0.4 

15 

15.00 

0.0 

0.0 

2.1 

43 

21.17 

0.0 

0.0 

0.5 

16 

15.17 

0.0 

0.1 

0.0 

44 

21.29 

0.0 

0.1 

0.3 

17 

15.56 

0.1 

0.1 

1.3 

45 

21.41 

0.0 

0.0 

0.6 

18 

16.26 

0.0 

0.0 

0.5 

46 

21.45 

0.0 

0.0 

0.0 

19 

16.60 

0.0 

0.0 

2.2 

47 

21.55 

0.0 

0.2 

0.0 

20 

16.75 

0.0 

0.0 

0.0 

48 

21.60 

0.0 

0.0 

0.3 

21 

17.11 

0.1 

0.1 

0.3 

49 

21.63 

0.2 

0.0 

0.1 

22 

17.26 

0.0 

0.0 

0.3 

50 

21.69 

0.0 

0.0 

0.8 

23 

17.38 

0.1 

0.0 

0.3 

51 

21.89 

0.0 

0.0 

0.1 

24 

17.77 

0.1 

0.0 

3.4 

52 

21.96 

0.0 

0.0 

0.0 

25 

17.95 

0.0 

0.1 

0.1 

53 

22.17 

0.0 

0.1 

0.5 

26 

18.27 

0.0 

0.0 

0.1 

54 

22.18 

0.3 

0.0 

0.0 

27 

18.51 

0.0 

0.0 

0.3 

55 

22.33 

0.0 

0.0 

0.1 

28 

18.68 

0.0 

0.0 

0.5 

56 

23.09 

0.0 

0.2 

5.5 
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Table 2. WP-02AVP-04 Combined Cargo Element Modal Survey Target Modes 
Pre-test Analysis Sets Before and After Considering Mode Participation 


Only Considering MEM 

After Considering Mode Participation 

Mode 

Freq. 

X„MEM 

Y^MEM 

ZoMEM 

Mode 

Freq. 

X„MEM 

Y„MEM 

Z„MEM 


(Hz) 

(%) 

(%) 

(%) 


(Hz) 

(%) 

(%) 

(%) 


6.66 


87.5 



6.66 


87.5 



8.30 

49.8 


21.4 


8.30 

49.8 


21.4 


9.40 

46.4 


26.7 


9.40 

46.4 


26.7 


10.60 


6.0 



10.60 


6.0 



13.08 



2.1 


13.08 



2.1 

12 

14.26 



4.7 

12 

14.26 



4.7 

13 

14.33 



8.6 

13 

14.33 



8.6 

14 

14.51 



3.7 

14 

14.51 



3.7 

15 

15.00 



2.1 

15 

15.00 



2.1 

19 

16.60 



2.2 

19 

16.60 



2.2 

24 

17.77 



3.4 

24 

17.77 



3.4 






28 

18.68 

0.0 

0.0 

0.5 






29 

19.09 

0.0 

0.0 

0.5 






30 

19.16 

0.0 

0.0 

1.4 

56 

23.09 



5.5 

56 

23.09 



5.5 

To 

tal 

97.2 

94.1 

80.5 

Total 

97.2 

94.1 

82.9 


Note; MEM percentages not shown are less than 2.0%; however, total MEM percentages include these values. 
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Figure 1.— The Space Station WP-(^/WP-04 Combined Cargo Element Space 
Shuttle payload. 


opttoniil 



Figure 2. — Mode participation analyses are performed as modal 
transient analysis data recovery restarts. 



I stop I 


Figure 3. — Operational flow for mode participation DMAP sequence. 
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Figure 4. — ^Finite element representation of the WP-02/WP-04 Combined 
Cargo Element for dynamic loads analyses. 


O 0 © 



© Collector 

Superelememt 



© Collector 

Superelement 

Figure 5. — Multi-level superelement tree for the WP-02/WP-04 
Combined Cargo Element dynamic loads analysis model. 




Figure 6. — WP-02/WP-04 Combined Cargo Element modes 28, 29, and 30 were 
identified as candidate modal survey target modes based on mode participation 
analysis. 
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DMAP Sequence for 
Mode Participation Analysis 
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■ -0> </>«■ </> </><*<«■</><*</>«■ -U> </> M </> -tA </> -W </> -W -tft -w </></> 


SUBDMAP MODEPART $ 

$ CALCULATE MODE PARTICIPATION 


(MODEPART) 


MODE PARTICIPATION FRACTIONS 
MSC/NASTRAN VERSION 67 

THIS DMAP CALCULATES MODE PARTICPATION FRACTIONS (MPF) FOR THE 
GENERALIZED (Q-SET) D.O.F. CONTRIBUTIONS OF A COMPONENT RESPONSE. 

THE MPF FOR A PARTICULAR Q-SET D.O.F. IS DEFINED AS THE PARTICULAR 
Q-SET D.O.F. CONTRIBUTION DIVIDED BY THE TOTAL Q-SET CONTRIBUTION TO 
THE TOTAL RESPONSE. THE RESPONSES ARE RECOVERED VIA DATA RECOVERY 
MATRICES (DRM) FORMED BY THE RUNS DRMl AND DRM2 . THE Q-SET D.O.F. 
BELONG TO THE COMPONENT DEFINED AS THE "BASE" SUPERELEMENT IN THE 
DRM RUNS. 


REQUIREMENTS TO USE THIS DMAP SEQUENCE - 


FILE MANAGEMENT SECTION - 

THE FOLLOWING DATA BASE FILES MUST BE AVAILABLE - 

- RUN DATA BASE 

- SYSTEM MODES DATA BASE 

- COMPONENT DATA BASE (IF SEPARATE FROM SYSTEM DATA BASE) 

- QTRAN DATA BASE 

- FORCING FUNCTION DATA BASE 

- FFMAX DATA BASE 

- UNTLD DATA BASE 

THE FOLLOWING FORTRAN FILES MUST BE AVAILABLE - 

- DRM FILE (CREATED VIA DRMl, DRM2, AND DRMEXT) 

NOTE! THE FILE MAY NOT BE NEEDED IF A RESTART IS REQUESTED 


EXECUTIVE CONTROL DECK - 

THIS DMAP SEQUENCE (DO NOT USE A RIGID FORMAT) 
DIAG 8, 14, AND 20 RECOMMENDED 


CASE CONTROL DECK - 

NO REQUIREMENTS - STANDARD TITLING (TITLE, SUBTITLE, LABEL) 


BULK DATA DECK - 


REQUIRED PARAMATER CARDS - 


PARAM, CSElU 

^r«rM-innT 'T'ktTT'i.lm Tr\ tPUC rtDM fPUTG 

— iTiii i>\jtrctr\C/Uc,L'iCii'< 1. xl/ rvxi\ x4»j-i 4..**^ 

SHOULD BE THE "BASE SUPERELEMENT' ID FROM THE 
RUNS (DRMl AND DRM2) 

(NO DEFAULT - REQUIRED) 

ID 

DRM 

PARAM, DLOAD 

- DEFINE THE LOAD CASE ID 
(NO DEFAULT - REQUIRED) 


PARAM, FFSEID 

- DEFINE THE FORCED SUPERELEMENT ID 
(NO DEFAULT - REQUIRED) 



IMPORTANT PARAMATER CARDS - ONE OF THE FOLLOWING MUST BE RESET 

PARAM.DRMACCE - ’-I" INDICATES DRM IS FOR ACCELERATIONS ONLY 
(DEFAULT = 0 - DRM NOT ONLY ACCE) 

PARAM,DRMALL - "-I" INDICATES DRM IS FOR ACCE AND DISP 
(DEFAULT = 0 - DRM NOT FOR ACCE AND DISP) 


PARAM.DRMDISP 


-1* INDICATES DRM IS FOR DISPLACEMENTS ONLY 


• </></></></></> </v </></></> </v ^ ^ ^ 


(DEFAULT = 0 - DRM NOT ONLY DISP) 


OPTIONAL PARAMATER CARDS - 

PARAM,DRMID - DRM ID USED WHEN STORING DRM ON DATA BASE 
(DEFAULT = 0) 

PARAM,DYNFCTR - FACTOR APPLIED TO THE ELASTIC PORTION OF THE 
DYNAMIC RESPONSE 
(DEFAULT =1.0) 

PARAM, FACTOR - "-1” INDICATES MODAL RESPONSES AND APPLIED 'LOADS 
ARE NOT BEING FACTORED (FASTER RUN) 

(DEFAULT = 0 - FACTOR EVEN THOUGH MAY BE 1.0) 

PARAM, FUNITID - FORTRAN UNIT ID OF DRM FILE 
(DEFAULT = 14) 

PARAM, NOMODACC - ”-1“ INDICATES USE MODE DISPLACEMENT FOR SOLVING 
BOUNDARY SOLUTION 

(DEFAULT = 0 - USE MODE ACCELERATION) 

PARAM, NORBDISP - "-1" INDICATES RIGID-BODY MODAL DISPLACEMENTS ARE 
SET TO ZERO 

(DEFAULT = 0 - KEEP MODAL RB DISPLACEMENTS) 

PARAM, RESPONSE - DEFINE THE RESPONSE NO. OF THE DRM FOR CALCULATING 
MODE PARTICIPATION FRACTIONS. 

(DEFAULT = -1 - MPF'S FOR ALL RESPONSES IN DRM) 

PARAM, RESTART - •-1" INDICATES RESTART RUN WHERE BOUNDARY SOLUTION 
AND DRM ARE FETCHED FROM THE DATA BASE. FRACTIONS 
ARE CALCULATED FOR A NEW RESPONSE AND/OR NEW TIME 
(DEFAULT = 0 - PERFORM ENTIRE RUN) 

PARAM, SOLTIME - DEFINE THE COLUMN NO. OF THE BOUNDARY SOLUTION 
CORRESPONDING TO TIME OF INTEREST. 

(DEFAULT = -1 - AUTOMATIC SEARCH COL OF ABS MAX) 

PARAM, STATFCTR - FACTOR APPLIED TO THE STATIC PORTION OF THE 
RESPONSE 
(DEFAULT =1.0) 

PARAM, TSETRESP - IF +1, PORTIONS OF TOTAL RESPONSE DUE TO T-SET DOF 
ACCELERATIONS AND DISPLACEMENTS PRINTED 
(DEFAULT = -1 - DO NOT PRINT T-SET RESPONSES) 


EXAMPLE NASTRAN DECK - 

ID MODE, PARTICIPATION 

ASSIGN MASTER = 'MODEPART. MASTER' 
ASSIGN DBALL = 'MODEPART. DBALL' 
ASSIGN USRSOU = ' MODEPART. USRSOU' 
ASSIGN USROBJ = ' MODEPART. USROBJ' 
INIT DBALL, LOGICAL= (DBALL (100000) ) 
$ 

ASSIGN FF=TRLG. MASTER $ 

DBLOCATE LOGI=FF $ 

$ 

ASSIGN FFX= FFMAX. MASTER $ 

DBLOCATE LOGI=FFX $ 

$ 

ASSIGN QTR=QTRAN. MASTER $ 

DBLOCATE LOGI=QTR $ 


$ ASSIGN UNIT=UNTLD. MASTER $ 

$ DBLOCATE LOGI=UNIT $ 

$ $ 

$ ASSIGN SYS=SySMODES. MASTER $ 

$ DBLOCATE LOGI=SYS $ 

$ $ 

$ ASSIGN INPUTT4=DRM.OPT4, UNIT=11, FORM-FORMATTED $ 

$ $ 

$ ACQUIRE NDDLOLD 

$ S 

$ TIME 60 $ 

$ DIAG 8,14,20 $ 

$ SOL MODEPART $ 

$ COMPILE MODEPART, NOLIST, NOREF 

$ INCLUDE MODEPART. V67 

$ CEND 

$ TITLE = MODE PARTICIPATION FRACTIONS 

$ SUBTITLE = ALL DRM RESPONSES (Tj-ME OF ABSOLUTE MAX) 

$ LABEL = LOAD CASE 2165 

$ $ 

$ BEGIN BULK 

$ $ 

$ $ PARAMATERS 

4 $ 

$ PARAM,CSEID,100 $ DEFINE DRM SUPERELEMENT ID 

$ PARAM,DLOAD,2165 $ DEFINE DLOAD ID 

$ PARAM,DRMALL,-1 $ DRM CONTAINS ACCE AND DISP 

$ PARAM,DRMID,101 $ DEFINE THE DRM ID (FOR DBSTORE) 

5 PARAM,DYNFCTR,1.50 $ DEFINE DYNAMIC FACTOR 

$ PARAM,FFSEID, 99 $ DEFINE FORCED SUPERELEMENT ID 

$ PARAM,STATFCTR=1.0 $ DEFINE STATIC FACTOR 

$ $ 

$ ENDDATA 

$ 

$ 

$ HISTORY DOCUMENTATION - 

$ 

$ VERSION 1.0 A. BARNETT & O. IBRAHIM ll-MAY-92 
S - ORIGINAL VERSION 

$ - CALCULATE PARTICIPATION FRACTIONS FOR ANY RESPONSE 

$ USING DATA RECOVERY MATRICES 

$ VERSION 1.1 NADINE MACK 05-JUN-92 

$ - CONVERTED TO VERSION 67 

$ VERSION 1.2 A. BARNETT & O. IBRAHIM 18-JUN-92 

$ - CORRECT LOGIC WHEN TOTAL DRM IS FORMED FROM SEPARATE 

$ PIECES READ AS INPUTT4 . IT IS ASSUMED THAT ALL T- 

$ AND Q-SET D.O.F. ARE REPRESENTED IN FINAL MATRIX. 

$ VERSION 1.3 A. BARNETT & O. IBRAHIM 14-MAY-93 

$ - ADD LOGIC FOR AUTOMATIC SEARCH FOR TIME OF ABS MAX 

$ RESPONSE AND CORRESPONDING MPF CALCULATION 

$ - ADD LOGIC FOR AUTOMATIC LOOPING OVER ALL RESPONSES 

$ WITHIN DRM 

$ VERSION 1.4 A. BARNETT & O. IBRAHIM 17-MAY-93 

$ - ADD OPTIONAL PARAM FOR PRINTING T-SET RESPONSES 

$23456789012345678901234567890123456789012345678901234567890123456789012 

$ 1 2 3 4 5 6 7 

$ 

$ SET DEFAULT VALUES 
$ 


TYPE 

PARM, 

,I,Y,CONTINU=-l $ 

TYPE 

FARM, 

,I,Y,CSEID=-1 $ 

TYPE 

PARM, 

,I,Y,DLOAD=0 $ 

TYPE 

PARM, 

,I,Y,DRMACCE=^0 $ 

TYPE 

PARM, 

• T.Y.DRMALL=0 $ 


CONTINUE RESPONSES 

DRM SUPERELEMENT ID 

DLOAD ID 

DRM FOR ACCE 

DRM FOR ACCE AND DISP 


TYPE 

FARM, ,I,Y,DRMDISP=0 $ 

DRM FOR DISP 

TYPE 

PARM, ,I,Y,DRMID=0 $ 

DRM ID 

TYPE 

PARM, , RS, Y, DYNFCTR=1 . 0 $ 

DYNAMIC FACTOR 

TYPE 

PARM, ,I,Y,EXTTIME=1 $ 

TIME TO EXTRACT 

TYPE 

PARM, ,I,Y,FACTOR=0 $ 

PERFORM FACTORING 

TYPE 

PARM, ,I,Y,FFSEID=-1 $ 

FORCED SUPERELEt-IENT ID 

TYPE 

PARM, ,I,Y,FUNITID=14 $ 

FORTRAN UNIT ID 

TYPE 

PARM , , I , Y , NOMODACC= 0 $ 

USE MODE ACCELERATION 

TYPE 

PARM, ,I,Y,NORBDISP=0 $ 

KEEP RB MODAL DISP 

TYPE 

PARM, , I , Y , RESPONSE= - 1 $ 

RESPONSE NO. FROM DRM 

TYPE 

PARM, , I , Y , RESTART=0 $ 

RESTART RUN 

TYPE 

PARM, , I,Y,SOLTIME=-l $ 

SOLUTION TIME 

TYPE 

PARM, , I , y , STARTCNT=0 $ 

COL TO START EXT SEARCH 

TYPE 

PARM, , I , Y, STARTRES=0 $ 

COL TO START RESPONSES 

TYPE 

PARM, , RS, Y, STATFCTR=1 . 0 $ 

STATIC FACTOR 

ITYPE 

PARM, , I , Y , TSETRESP=-1 $ 

PRINT T-SET RESrONSES 

[type 

PARM,,I,N,CARDNO=0,PFILE=0 $ 

CARD AND PLOT COUNTERS 

$ PATH QUALIFIERS 
$ 


$ 

TYPE 

PARM,NDDL,I,Y,PEID $ 


TYPE 

PARM,NDDL,I,N,SEID $ 


§ 

$ ALLOW 
$ 

OVERWRITES 


$ 

FILE 

DRMQX=OVRWRT $ 


FILE 

DRMTAX=OVRWRT $ 


FILE 

DRMTDX=OVRWRT $ 


FILE 

DRMTX=OVRWRT $ 



FILE DRMX=OVRWRT $ 

FILE NOQUNIT=OVRWRT $ 

FILE NOQUNITT=OVRWRT $ 

FILE QSETMPF=OVRWRT $ 

FILE QSETMPFT=OVRWRT $ 

FILE QSETRESP=0VRWRT $ 

FILE QSETRSPX=OVRWRT $ 

FILE QSETTERM=0VRWRT $ 

FILE QSTTERMT=OVRWRT $ 

FILE SOLN=OVRWRT $ 

FILE TOTRESP=OVRWRT $ 

FILE TOTSOLX=OVRWRT $ 

FILE TSETRESP=OVRWRT $ 

FILE TSETRSPA=OVRWRT $ 

FILE TSETRSPD=OVRWRT $ 

FILE UQDDX=OVRWRT $ 

FILE UQDDXT=OVRWRT $ 

FILE UTDDXsOVRWRT $ 

FILE UTX=OVRWRT $ 

FILE VSOLNMX=OVRWRT $ 

$ 

$ ASSIGN DATA BLOCK DEFINITIONS 
$ - 


$ 

TYPE 

DB,EQEXINS 

$ 

PATH=MP 

LOCATION=DBUP 

TYPE 

DB,GPLS 

$ 

PATH=MP 

LOCATION=DBUP 

TYPE 

DB,SILS 

$ 

PATH=MP 

LOCATION=DBUP 

TYPE 

DB,USET 

$ 

PATH=MP 

LOCATION=DBDN 


$ 

$ DMAP DATA BLOCK AND PARAMETER ALLOCATION TO DATA BASE 

$ THE PATHS ARE DEFINED AS FOLLOWS: 

$ PATH MP UZR, MODEL, PEI D $ 


TYPE, PARM.NDDL, CHARS, N,IFPX $ LOCATIONS 

type, PARM.NDDL.I.Y.M0DEL,NL99,SEIDP,S0LID.S0LP,UZR $ QUALIFIERS 


* 

TYPE DB CASECC,PCDB,XYCDB,POSTCDB, FORCE, BULK, EPT, MPT, EDT.DIT, 

DYNAMICS , GEOMl , GBOM2 , GEOM3 , GEOM4 , MATPOOL , AXIC, PVT, DMI , 

DTI, DTINDX, EDOM $ 

TYPE PARM, , CHARS, Y,IFP=' DEFAULT' 

$ 

IFPl /CASECC, PCDB, XYCDB, POSTCDB, FORCE/S, N, NOGOIFPl $ 

XSORT FORCE, /IBULK/S,N,NOGOXSRT $ FOR NODATA ETC. OPTION 
IFP IBULK/ 

GEOMIQ, lEPT, IMPT, lEDT, IDIT, IDYNAMIC, IGEOM2 , IGEOM3 , IGEOM4, , 

FP6SCR, IMATPOOL,AXIC, PVT, DMI, DMINDX, DTI, DTINDX, , 

EDOM, DEQATN , DEQINDX/ 

S,N,NOGOIFP $ 

IF (NOGOIFPl OR NOGOXSRT OR NOGOIFP) THEN $ 

MESSAGE //'ERROR IN USER INPUT - RUN TERMINATED'/ S 
EXIT S 
ENDIF $ 

$ 

IF (IFP o'DEFAULT') IFPX=IFP $ 

« 

SETVAL //NP/S, Y,MODEL/NP/S,Y,NL99/NP/S,Y,SEIDP $ QUALIFIERS 
SETVAL //NP/S, Y,SOLID/NP/S,Y,SOLP/NP/S,Y,UZR $ QUALIFIERS 

PARAMR //'COMPLEX' //STATFCTR/0. 0/S, N,STATFCTC $ COMPLEX FOR ADD MOD 

PARAMR //'COMPLEX'//DYNFCTR/0.0/S,N,DYNFCTC $ COMPLEX FOR ADD MODULE 

$ 

$ DATABASE IS VALID FROM 'THIS POINT ON 

5 - 

$ 

PUTSYS(1,125) $ HONOR RESTART/NOKEEP OPTION FROM NOW ON. 

S 

COND NODBPRT, DBDICT $ PRINT DATA BASE DICTIONARY 

DBDIR //V,Y,DBDRPRJ=0/V,Y,DBDRVER=0/V,Y,DBDROPT=63 $ 

PARAM //'NOP'/S,Y,DBDICT=-l $ SET DEFAULT 

LABEL NODBPRT $ 

g 

$ CHECK THAT THE SUPERELEMENT ID IS SPECIFIED 
$ - 

PARAM //'LE'/V,N,BADSEID/CSEID/0 $ -1 IF CSEID . LE . 0 

IF (BADSEID=-1) THEN $ 

MESSAGE //'FATAL ERROR - CSEID HAS NOT BEEN SET' / $ 

ENDIF $ 

$ 

$ CHECK THAT THE DLOAD SET ID IS SPECIFIED 

S 

$ 

PARAM //'LE'/S,N,BADDLOAD/DLOAD/0 $ -1 IF DLOAD .LE. 0 

IF (BADDLOAD=-l) THEN $ 

MESSAGE //'FATAL ERROR - DLOAD HAS NOT BEEN SET'/ $ 

ENDIF S 

$ 

$ CHECK THAT THE FORCED SUPERELEMENT IS SPECIFIED 

$ 

PARAM //'LE'/S,N,BADFFID/FFSEID/0 $ -1 IF FFSEID .LE. 0 

IF (BADFFID=-1) THEN $ 

MESSAGE //'FATAL ERROR - FFSEID HAS NOT BEEN SET'/ $ 

ENDIF $ 

$ 

$ TERMINATE THE RUN IF ANY ERRORS ENCOUNTERED 

5 

4 

COND ERROR, BADDLOAD $ ERROR IF BAD DLOAD 

COND ERROR, BADFFID $ ERROR IF BAD FFSEID 

COND ERROR, BADSEID 9 ERROR IF BAD SEID 


FETCH SUPERELEMENT COMPONENT DATA 


$ 

DBVIEW CSILS = SILS (WHERE SEID=CSEID) $ 

DBVIEW CGPLS = GPLS (WHERE SEID=CSEID) $ 

DBVIEW CUSET = USET (WHERE SEID=CSEID) $ 

DBVIEW CEQEXINS = EQEXINS (WHERE SEID=CSEID) $ 

PARAML CUSET/ / ' USET' / /////'A'/S,N, NOASET/ 

'T'/S,N,NOTSET/ 

'Q'/S,N,NOQSET $ NO. OF DOF 

$ 

$ IF THIS IS A RESTART, GO TO MPF CALCULATION 

$ 

$ 

COND CALCMPF, RESTART $ 

$ 

$ EXTRACT THE RIGID-BODY AND ELASTIC MODAL ACCELERATIONS 

$ 

$ 

CALL DBFETCH /UHVF, , , , /DLOAD/0 /O /O /O $ FETCH MODAL RESPONSES 

PARAML UHVF//'PRESENCE'////S,N,NOUHVF $ -1 IF NO UHVF 

IF (NOUHVF=-l) THEN $ 

MESSAGE //'FATAL ERROR - UHVF DATA BLOCK REQUIRED BUT NOT FOUND'/ $ 
ENDIF $ 

COND ERROR, NOUHVF $ ERROR IF NO UHVF 

PARAML UHVF//'TRAILER'/1/S,N,NOCOLU $ NO. OF COL IN UHVF 

PARAML UHVF/Z'TRAILER' /2/S, N, NOMODES $ NO. OF MODES IN UHVF 

MATGEN , /VUXA/4/1/NOCOLU//1/3/3/ $ ONES FOR ACCELERATIONS 

MATCEN , /VHRBEL/6/NOMODES/6/NOMODES $ H = RB / EL 

PARTN UHVF,VUXA,VHRBEL/,,QDDRB1,QDDEL1/1 $ ACCE OF RB AND EL 


IF THIS IS A RESTART, GO TO MPF CALCULATION 


CALCMPF, RESTART $ 


COND ERROR, NOUHVF $ ERROR IF NO UHVF 

PARAML UHVF//'TRAILER'/1/S,N,NOCOLU $ NO. OF COL IN UHVF 

PARAML UHVF//'TRAILER' /2/S, N, NOMODES $ NO. OF MODES IN UHVF 

MATGEN , /VUXA/4/1/NOCOLU//1/3/3/ $ ONES FOR ACCELERATIONS 

MATCEN , /VHRBEL/6/NOMODES/6/NOMODES $ H = RB / EL 

PARTN UHVF,VUXA,VHRBEL/,,QDDRB1,QDDEL1/1 $ ACCE OF RB AND EL 

g 

$ APPLY FACTORS TO THE RIGID-BODY AND ELASTIC MODAL ACCELERATIONS 


PARAM, STATFCTR APPLIED TO RIGID-BODY MODE CONTRIBUTIONS 
PARAM, DYNFCTR APPLIED TO ELASTIC MODE CONTRIBUTIONS 


QDDRB1,QDDRB/FACTOR $ 
QDDELl,QDDEL/FACTOR $ 
DNFACTRl, FACTOR $ 

QDDRBl, /QDDRB/STATFCTC $ 
QDDELl, /QDDEL/DYNFCTC $ 
DNFACTRl $ 


EQUIV IF NO FACTORING 
EQUIV IF NO FACTORING 
SKIP IF NO FACTORING 
MULT BY STATIC FACTOR 
MULT BY DYNAMIC FACTOR 
END OF FACTORING 


LABEL DNFACTRl S 

$ 

$ FORM COMPONENT T-SET / Q-SET PARTITIONING VECTOR 

g 

g 

VEC CUSET/CVATQ/'A'/'COMP'/'Q' $ ONES FOR Q-SET DOF 

MATPRN CVATQ// $ PRINT 

MATCPR CGPLS, CUSET, CSILS, CVATQ//'H'/'A' $ PRINT 

$ 

$ PARTITION THE COMPONENT SYSTEM MODES 

g 

$ 

$ THEY ARE FIRST PARTITIONED INTO RIGID-BODY AND ELASTIC SETS; THEN 
$ ACCORDING TO T-SET AND Q-SET DOF MEMBERSHIP 

g 

CALL DBFETCH /HULVS, , , , /O/CSEID/1/0/0 $ A-SET MODE SHAPES 

PARAML HULVS//' PRESENCE' ////S,N,BADHULVS $ -1 IF NO ULVS 

IF (BADHULVS=-1) THEN $ * 

MESSAGE //'FATAL ERROR - ULVS FOR MODES REQUIRED BUT NOT FOUND'/ $ 
ENDIF ; 


THEY ARE FIRST PARTITIONED INTO RIGID-BODY AND ELASTIC SETS; THEN 
ACCORDING TO T-SET AND Q-SET DOF MEMBERSHIP 

DBFETCH /HULVS, ,,, /O/CSEID/1/0/0 $ A-SET MODE SHAPES 

ML HULVS//'PRESENCE'////S,N,BADHULVS $ -1 IF NO ULVS 



COND ERROR, BADHULVS $ 

PARTN HULVS , VHRBEL, /HULVSRB, , HULVSEL, /I 

PARTN HULVSRB, , CVATQ/PHITR, PHIQR, , /I $ 
PARTN HULVSEL, ,CVATQ/PHITE, PHIQE, , /I $ 

$ 

$ CALCULATE COMPONENT T-SET, Q-SET, AND A- 

$ 

$ 

mpyad phitr,qddrb,/utddrb $ 

mpyad phite,qddel,/utddel $ 

add utddrb,utddel/utdd $ 

MPYAD PHIQR, QDDRB,/UQDDRB $ 

MPYAD PHIQE, QDDEL,/UQDDEL $ 

add UQDDRB,UQDDEL/UQDD $ 

MERGE UTDD, UQDD, , , , CVATQ/UADD/1 $ 

$ 

$ CALCULATE COMPONENT T-SET DISPLACEMENTS 


ERROR IF NO ULVS 
$ COL PARTITION 
ROW PARTITION 
ROW PARTITION 

SET ACCELERATIONS 


UTDDRB = PHITR * QDDRB 
UTDDEL = PHITE * QDDEL 
UTDD = UTDDRB + UTDDEL 
UQDDRB = PHIQR * QDDRB 
UQDDEL = PHIQE * QDDEL 
UQDD = UQDDRB + UQDDEL 
ROW MERGE 

VIA MODE ACCELERATION 


MATGEN , /VUXD/4/1/NOCOLU//1/3/1/ $ ONES FOR DISPLACEMENTS 

PARTN UHVF, VUXD, VHRBEL/, ,QRB1,QEL1 /I $ DISP OF RB AND EL 

COND MDDISP,NOMODACC $ SKIP IF MODE DlSP 

$ 

$ EXTRACT THE NON-ZERO ROWS FROM THE APPLIED LOADS MATRIX 

5 - 

$ 

CALL DBFETCH /PDFX, , , , /DLOAD/FFSEID/0/0/0 $ APPLIED LOADS 

PARAML PDFX//'PRESENCE'////S,N,BADPDF $ -1 IF NO PDF 

IF (BADPDF=-1) THEN $ , . 

MESSAGE //'FATAL ERROR - PDF DATA BLOCK REQUIRED BUT NOT FOUND'/ $ 
ENDIF $ 

COND ERROR, BADPDF $ ERROR IF NO PDF 

CALL DBFETCH /PGMAXALL, , , , /O/FFSEID/0/0/0 $ MAX. DYNAMIC LOADS 

PARAML PGMAXALL/ /' PRESENCE '////S,N,BADPGMAX $ -1 IF NO PGMAXALL 

IF (BADPGMAX=-1) THEN $ 

MESSAGE //'FATAL ERROR - PGMAXALL REQUIRED BUT NOT FOUND'/ $ 

ENDIF $ 

COND ERROR, BADPGMAX $ ERROR IF NO PGMAXALL 

PARTN PDFX, , PGMAXALL/, PDFNZl, , /I $ KEEP NON- ZERO ROWS 

$ 

S APPLY FACTOR TO THE STATIC PORTION OF THE RESPONSES 

5 - - - 

$ 

$ PARAM,STATFCTR APPLIED TO UNIT IXJAD CONTRIBUTIONS 

EQUIV PDFNZl, PDFNZ/FACTOR $ EQUIV IF NO FACTORING 

COND DNFACTR2 , FACTOR $ SKIP IF NO FACTORING 

ADD PDFNZl, /PDFNZ/STATFCTC $ MULT BY STATIC FACTOR 

LABEL DNFACTR2 $ END OF FACTORING 

$ 

$ FORM THE INVERSE OF THE EIGENVALUE (OMEGA-SQ) MATRIX 

$ 

$ 

CALL DBFETCH /KHH, MHH, , , /DLOAD/0/0/0/0 $ SYSTEM DATA 

ADD MHH,KHH/OMSQINVl///2 $ MATRIX "DIVIDE” 

PARTN OMSQINVl, VHRBEL, /,, ,OMSQINV/ $ ELIM. RB MODES 

$ 

$ PARTITION 'THE UNIT LOAD STATIC DEFLECTIONS TO T-SET DOF 

5 

$ 

CALL DBFETCH /SULVS, , , , /O/CSEID/0/0/0 $ A-SET DEFLECTIONS 

PARAML SULVS/Z'PRESENCE' ////S,N,BADSULVS $ -1 IF NO ULVS 

IF (BADSULVS = -1) THEN $ 

MESSAGE //'FATAL ERROR - SULVS DATA BLOCK REQUIRED BUT NOT FOUND'/ $ 
ENDIF $ 

COND ERROR, BADSULVS S ERROR IF NO ULVS 


PARTN SULVS, ,CVATQ/SULVST, , ,/l $ ROW PARTITION 

$ 

$ CALCULATE T-SET DISPLACEMENTS DUE TO ELASTIC RESPONSE 

$ 

$ 

SMPYAD PHITE, OMSQINV, QDDEL, , ,/UTEL/3 /-I $ -l*PHITE*OMSQINV*QDDEL 

MPYAD SULVST, PDFNZ, /UTPDF $ UTPDF = SULVST * PDFNZ 

ADD UTEL.UTPDF/UTl S UTl = UTEL + UTPDF 

JUMP DNUTl $ goto TO RB DISP 

$ 

$ CALCULATE COMPONENT T-SET DISPLACEMENTS VIA MODE DISPLACEMENT 


$ 

LABEL MDDISP $ 

$ 

$ APPLY FAC'TOR TO THE ELASTIC DYNAMIC PORTION OF THE RESPONSES 


EQUIV IF NO F.ACTORING 
SKIP IF NO FACTORING 
MULT BY DYNAMIC FACTOR 
END OF FACTORING 


$ 

$ PARAM,DYNFCTR APPLIED TO ELASTIC MODE CONTRIBUTIONS 

$ 

EQUIV QELl,QEL/FACTOR $ EQUIV IF NO F.AC'TORING 

COND DNFACTR3 , FAC'TOR $ SKIP IF NO FACTORING 

ADD QELl, /QEL/DYNFC'TC $ MULT BY DYNAMIC FACTOR 

LABEL DNFACTR3 $ END OF FACTORING 

$ 

$ CALCULA'TE T-SET DISPLACEMENTS DUE TO ELASTIC RESPONSE 

$ - 

$ 

MPYAD PHITE, QEL, /UTl $ UTl = PHITE * QEL 

$ 

$ INCLUDE RIGID-BODY CONTRIBUTION TO T-SET DISPLACEMENTS (IF REQUIRED) 


QELl,QEL/FACTOR $ 
DNFACTR3 , FAC'TOR $ 
QELl, /QEL/DYNFC'TC $ 
DNFACTR3 $ 


DNUTl $ 

UTl,UT/NORBDISP $ 
DNSOLVE,NORBDISP $ 


EQUIV IF NO RB DISP 
SKIP IF NO RB DISP 


$ APPLY FACTOR TO THE RIGID-BODY DYNAMIC PORTION OF THE RESPONSES 

s - — 

$ 

$ PARAM,STATFCTR APPLIED TO RIGID-BODY MODE CONTRIBUTIONS 

$ 

EQUIV QRBl,QRB/FACTOR $ EQUIV IF NO FACTORING 

COND DNFACTR 4, FACTOR $ SKIP IF NO FACTORING 

ADD QRBl, /QRB/DYNFCTC $ MULT BY DYNAMIC FACTOR 

LABEL DNFACTR4 $ END OF FACTORING 

$ 

$ CALCULATE T-SET DISPLACEMENTS DUE TO RIGID-BODY RESPONSE 

$ 

$ 

MPYAD PHITR, QRB,/UT2 $ UT2 = PHITR * QRB 

ADD UT1,UT2/UT $ UT = UTl + UT2 

$ 

S STORE ACCELERATIONS AND DISPLACEMENTS ON DATA BASE 

$ 

$ 

LABEL DNSOLVE $ 

CALL DBSTORE UADD, UTDD, UQDD, ,/ /DLOAD/CSEID/ ' DBALL' /O $ STORE ACCELS 

CALL DBSTORE UT, ,,, //DLOAD/CSEID/ ' DBALL' /O $ STORE DISPLACEMENTS 

$ 

$ MERGE THE COMPONENT BOUNDARY RESPONSES 

$ 

s 

PARAM //'ADD'/S,N,NOAT/NOASET/NOTSET $ NOAT = NOASET + NOTSET 

MATGEN , /VAT/6/NOAT/NOASET/NOTSET $ ZEROES FOR A-SET 

MERGE UADD.UT VAT/UADDUT/1 S ROW MERGE 


PARAM, STATFCTR APPLIED TO RIGID-BODY MODE CONTRIBUTIONS 


QRBl,QRB/FACTOR $ 
DNFACTR 4, FACTOR $ 
QRBl, /QRB/DYNFCTC $ 
DNFACTR4 $ 


EQUIV IF NO FACTORING 
SKIP IF NO FACTORING 
MULT BY DYNAMIC FACTOR 
END OF FACTORING 


DNSOLVE $ 

DBSTORE UADD, UTDD, UQDD, ,/ /DLOAD/CSEID/ 'DBALL' /O $ STORE ACCELS 
DBSTORE UT,,,, //DLOAD/CSEID/ 'DBALL' /O $ STORE DISPLACEMENTS 


$ 

$ FETCH THE DRM 

$ 

$ 

$ FOR ACCELERATIONS, THE DRM EQUATION IS 

$ 

$ ACCE = (ACCELRM * UADD) 

* 

$ FOR DISPLACEMENTS (DISP. DEPENDENT RESPONSES), THE DRM EQUATION IS 

$ 

$ DISP = (DISPLRM2 * UADD) + (DISPLRMl » UT) 

$ 

$ FOR ACCELERATIONS AND DISPLACEMENTS, THE DRM MATRIX EQUATION IS 

$ 

$ ACCE ACCELRM s 0 ^ UADD 

$ DISP DISPLRM2 : DISPLRMl UT 

$ 

$ A MAXIMUM OF FIVE DRM'S WILL BE READ. THEY ARE 
$ 1. CONTROL MATRIX (ALWAYS READ) 

$ 2 . ACCELRM 

$ 3 . DISPLRMl 

$ 4 . DISPLRM2 

$ 5. CONTROL MATRIX (ALWAYS READ) 

e 

$ THE FORMAT OF THE DRM INPUTT4 STATEMENT LINE DEPENDS UPON WHAT IS 
$ BEING DATA RECOVERED. 

COND ACCEONLY,DRMACCE $ GOTO ONLY ACCE DRM 

COND DISPONLY,DRMDISP $ GOTO ONLY DISP DRM 

COND ACCEDISP, DRMALL S GOTO ALL DRM 

MESSAGE //'FATAL ERROR - DRM TYPE HAS NOT BEEN SPECIFIED VIA PARAM' / $ 
JUMP ERROR $ ERROR IF NOTHING RESET 

$ 

LABEL ACCEONLY $ ^ 

INPUTT4 /DUM1,ACCELRM,DUM2,,/3/FUNITID//0 $ BCD READ ONLY ACCELRM 

PARAML ACCELRM// 'TRAILER' /1/S, N,NOCOLAC $ NO. OF COL OF ACCELRM 

PARAM //'EQ'/S,N,NOCOLEQA/NOCOLAC/NOASET $ -1 IF NOCOLAC = NOASET 

EQUIV ACCELRM, ACCELRMX/NOCOLEQA $ ACCELRMX = ACCELRM 

COND DNEXPNDA, NOCOLEQA $ SKIP IF NO EXPAND 

MATGEN , /VDRMNULL/6/NOASET/NOCOLAC/NOASET $ ZEROS FOR ACCELRM 

MERGE ACCELRM, ,, ,VDRMNULL,/ACCELRMX/1 $ COL MERGE W/ NULLS 

LABEL DNEXPNDA $ 

MERGE ACCELRMX, ,, ,VAT, /DRM/1 $ COL MERGE W/ NULLS 

MERGE UADD, , , , ,VAT/TOTSOL/l $ ROW MERGE W/ NULLS 

JUMP DNDRM $ 

$ 

LABEL DISPONLY $ 

INPUTT4 /DUMl, DISPLRMl, DISPLRM2, DUM2, /4/FUNITID//0 $ BCD ONLY DISP'S 
PARAML DISPLRM2//'TRAILER'/1/S,N,N0C0LD2 $ NO. OF COL OF DISPLRM2 

PARAM //'EQ'/S,N, NOCOLEQA /NOCOLD2 /NOASET $ -1 IF NOCOLD2 = NOASET 

EQUIV DISPLRM2,DISPLRMX/NOCOLEQA $ DISPLRMX = DISPLRM2 

COND DNEXPNDD, NOCOLEQA $ SKIP IF NO EXPAND 

MATCEN , /VDRMNULL/ 6 /NOASET/NOCOLD2 /NOASET $ ZEROS FOR DISPLRM2 

MERGE DISPLRM2, , , ,VDRMNULL, /DISPLRMX/1 $ COL MERGE W/ NULLS 

LABEL DNEXPNDD $ 

MERGE DISPLRMX, , DISPLRMl, , VAT, /DRM/1 $ DRM = DISP2 : DISPl 

EQUIV UADDUT, TOTSOL/ALWAYS $ TOTSOL = UADDUT 

JUMP DNDRM $ 

$ 

LABEL ACCEDISP $ , 

INPU'IT4 /DUMl, ACCELRM, DISPLRMl, DISPLRM2,DUM2/5/FUNITID//0 $ BCD ALL 
PARAML ACCELRM//'TRAILER'/2/S,N,N0ACCE $ NO. ROWS ACCELRM 

PARAML DISPLRM1//'TRAILER'/2/S,N,N0DISP $ NO. ROWS DISPLRMl 

PARAM //'ADD'/S,N,NOAD/NOACCE/NODISP S NOAD = NOACCE + NODISP 

MATGEN . /VAD/6/NOAD/NOACCE/NODISP $ ZEROES FOR ACCE 


ACCELRM// 'TRAILER' /1/S, N, NOCOLAC $ NO. OF COL OF ACCELRM 

//'EQ'/S,N,NOCOLEQA/NOCOLAC/NOASET $ -1 IF NOCOLAC = NOASET 

ACCELRM, ACCELRMX/NOCOLEQA $ ACCELRMX = ACCELRM 

DISPLRM2,DISPLRMX/NOCOLEQA $ DISPLRMX = DISPLRM2 

DNEXPND, NOCOLEQA $ SKIP IF NO EXPAND 

, /VDRMNULL/ 6 /NOASET/NOCOLAC /NOASET $ ONES FOR NULLS 
ACCELRM,,,, VDRMNULL, /ACCELRMX/1 $ COL MERGE W/ NULLS 

D1SPLRM2, ,, .VDRMNULL, /DISPLRMX/1 $ COL MERGE W/ NULLS 

DNEXPND $ 

ACCELRMX, DISPLRMX, , DISPLRMl , VAT, VAD/DRM/1 $ NON-SYM MERGE 
UADDUT, TOTSOL/AIWAYS $ TOTSOL = UADDUT 


DNDRM $ 

DBSTORE DRM, TOTSOL, 
SKPFTCH $ 


, //DLOAD/DRMID/'DBALL'/O $ STORE IN DBASE 
SKIP FETCH IF COLD RUN 


$ 

$ BEGIN CALCULATING MODE PARTICIPATION FRACTIONS 

$ 

$ 

LABEL CALCMPF $ * 

CALL DBFETCH /DRM, TOTSOL, ,, /DLOAD/DRMID/0/0/0 $ FETCH DRM & SOLN 

CALL DBFETCH /UQDD, UTDD, UT, , /DLOAD/CSEID/0/0/0 $ FETCH UQDD,UTDD,UT 

$ 

$ EXTRACT THE REQUIRED ROW FROM DRM 

$ 

$ 

LABEL SKPFTCH $ HERE IF COLD RUN 

IF (RESPONSE > 0) THEN $ SET FLAG FOR CONTINUE 

CONTINU = 0 $ 

ENDIF $ 

TRNSP DRM/DRMT $ TRANSPOSE DRM 

MATGEN , /VQDDX/6/NOAT/NOTSET/NOQSET/NOTSET $ ONES FOR Q-SET 

COND NOTSETl.TSETRESP $ JUMP IF NO T-SET RESP 

MATGEN , /VTDDX/6/NOAT/0/NOTSET/NOASET $ ONES FOR T-SET 

MATGEN , /VTX/6/NOAT/NOASET/NOTSET $ ONES FOR T-SET 

LABEL NOTSETl $ 

PARAML DRMT//'TRAILER'/l/S,N,NORESP $ NO OF RESPONSES 

PARAM //'ADD'/S,N,CNTRESP/STARTRES/1 $ BEGIN COUNTER 

LABEL TOPLOOPl $ TOP OF RESPONSE LOOP 

COND AUTGRESP, RESPONSE S GOTO ALL RESP SELECT 

JUMP GETRESP $ SKIP ALL RESP SELECT 

LABEL AUTORESP $ 

RESPONSE = CNTRESP $ SET RESPONSE NUMBER 

LABEL GETRESP $ 

MATMOD DRMT,,,,,/DRMTX, /1/RESPONSE $ EXTRACT COL FROM DRMT 

TRNSP DRMTX/DRMX $ ROW FROM DRM 

PARTN DRMX, VQDDX, /, ,DRMQX, /I $ COL PARTITION 

COND NOTSET2 , TSETRESP $ SKIP IF NO T-SET RESP 

PARTN DRMX.VTDDX, /,,DRMTAX, /I $ COL PARTITION 

PARTN DRMX.VTX, /,,DRMTDX, /I $ COL PARTITION 

LABEL NOTSET2 $ 

$ 

$ SET THE TIME INSTANT FOR MPF CALCULATION 


HERE IF COLD RUN 
SET FLAG FOR CONTINUE 


TRANSPOSE DRM 
ONES FOR Q-SET 
JUMP IF NO T-SET RESP 
ONES FOR T-SET 
ONES FOR T-SET 

NO OF RESPONSES 
BEGIN COUNTER 
TOP OF RESPONSE LOOP 
GOTO ALL RESP SELECT 
SKIP ALL RESP SELECT 

SET RESPONSE NUMBER 

EXTRACT COL FROM DRMT 

ROW FROM DRM 

COL partition 

SKIP IF NO T-SET RESP 

COL PARTITION 

COL PARTITION 


IF (SOLTIME > 0) THEN $ 

EXTTIME = SOLTIME $ 

JUMP GETSOLN $ 

ENDIF $ 

MPYAD DRMX, TOTSOL, /SOLN $ 

MATMOD SOLN, , , , , /VSOLNMX, /6 $ 

PARAML VSOLNMX/ /'DMI' /I /l/S,N,SOLNMX $ 

PARAM / / ' ADD ' /S , N , CNTTI ME /STARTCNT/ 1 $ 

LABEL TOPLOOP2 $ 

PARAML SOLN//'DMI'/CNTTIME/l/S,N,CMPARX $ 
PARAMR / / ' ABS ' /S , N , COMPARE /CMPARX $ 


DEFINE TIME TO EXTRACT 
SKIP AUTO TIME SELECT 

'CALCULATE SOLUTION 
ABSOLUTE MAX EXTREMA 
CONVERT VECTOR TO REAL 
BEGIN COUNTER 
TOP OF LOOP 
EXTRACT VALUE CNTTIME 
ABSOLUTE VALUE 


IF (COMPARE = SOLNMX) THEN $ 

EXTTIME = CNTTIME $ 

JUMP GETSOLN $ 

ENDIF $ 

CNTTIME = CNTTIME + 1 $ 

KEPT TOPLOOP2, 10000 $ 

$ 

$ EXTRACT THE REQUIRED COL FROM SOLUTIONS 

$ 

$ 

LABEL GETSOLN $ 

MATMOD TOTSOL, , , . » /TOTSOLX, /1/EXTTIME $ 
MATMOD UQDD, /UQDDX, /I /EXTTIME $ 

COND N0TSET3,TSETRESP $ 

MATMOD UTDD, , , , , /UTDDX, /1/EXTTIME $ 

MATMOD UT,,,,,/UTX, /1/EXTTIME $ 

LABEL NOTSET3 $ 


SET SOLUTION TIME 
JUMP TO GET REQ'D SOLN 

NEXT TIME STEP 


LABEL GETSOLN S 

MATMOD TOTSOL, /TOTSOLX, /1/EXTTIME $ EXTRA( 

MATMOD UQDD, ,,,, /UQDDX, / 1/EXTTIME $ EXTRA( 

COND N0TSET3,TSETRESP $ SKIP 

MATMOD UTDD, ,,,, /UTDDX, /1/EXTTIME $ EXTRA' 

MATMOD UT, ,,,,/UTX,/l/EXTTIME $ EXTRA' 

LABEL NOTSET3 $ 

$ 

$ CALCULATE THE TOTAL RESPONSE AND THE Q-SET RESPONSE 

$ 

$ 

MPYAD DRMX, TOTSOLX, /TOTRES.P $ TOTAL 

MPYAD DRMQX, UQDDX, /QSETRESP $ QSETR 

$ 

$ CALCULATE Q-SET MODE PARTICIPATION FRACTIONS 

$ 

$ 

TRNSP UQDDX/UQDDXT $ TRANS 

ADD DRMQX,UQDDXT/QSTTERMT///1 $ TERMW 

MATGEN , /NOQUNIT/6/NOQSET/0/NOQSET $ UNIT 

TRNSP NOQUNIT/NOQUNITT $ TRANS 

MPYAD QSETRESP, NOQUNITT,/QSETRSPX $ FULL 

ADD QSTTERMT,QSETRSPX/QSETMPFT///2 $ TERMW 

$ 

$ CALCULATE THE T-SET RESPONSE (IF REQUIRED) 

$ 

$ 

COND N0TSET4,TSETRESP $ SKIP 

MPYAD DRMTAX, UTDDX, /TSETRSPA $ T-SET 

MPYAD DRMTDX,UTX, /TSETRSPD $ T-SET 

ADD TSETRSPA, TSETRSPD/TSETRESP $ T-SET 

LABEL NOTSET4 $ 

$ 

$ OUTPUT RESULTS 

$ 


$ 

TRNSP QSTTERMT/QSETTERM $ 
TRNSP QSETMPFT/QSETMPF $ 
MATPRT TOTRESP// $ 

COND N0TSET5,TSETRESP $ 

MATPRT TSETRSPA// $ 

MATPRT TSETRSPD// $ 

MATPRT TSETRESP// $ 

LABEL NOTSET5 $ 

MATPRT QSETRESP// $ 

MATPRT QSETTERM/ / $ 

MATPRT QSETMPF// $ 

$ 

S NEXT RESPONSE (IF REQUIRED) 

j 

$ 

COND NEXTRESP, CONTINU $ 
JUMP FINIS $ 

LABEL NEXTRESP $ 

IF (CNTRESP = MORESP) THEN $ 

JUMP FINIS S 


EXTRACT COL FROM SOL 
EXTRACT COL FROM UQDD 
SKIP IF NO T-SET RESP 
EXTRACT COL FROM UTDD 
EXTRACT COL FROM UT 


TOTAL = DRMX * TOTSOLX 
QSETR = DRMQX * UQDDX 


TRANSPOSE UQDDX 
TERMWISE MULTIPLY 
UNIT VECTOR 
TRANSPOSE NOQUNIT 
FULL QSETRESP 
TERMWISE DIVIDE 


SKIP IF NO T-SET RESP 
T-SET ACCE RESP 
T-SET DISP RESP 
T-SET TOTAL RESP 


TRANSPOSE FOR OUTPUT 
TRANSPOSE FOR OUTPUT 
PRINT TOTAL RESPONSE 
SKIP IF NO T-SET RESP 
PRINT T-SET ACCE RESP 
PRINT T-SET DISP RESP 
PRINT T-SET TOTAL RESP 

PRINT Q-SET RESPONSE 
PRINT Q-SET TERMS 
PRINT Q-SET MPF'S 


TO END IF NOT ALL RESP 


TO END IP DONE 
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